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A simple model of spacetime foam, made by N Reissner-Nordstrom wormholes with a magnetic 
and electric charge in a semiclassical approximation, is taken under examination. The Casimir- 
like energy of the quantum fluctuation of such a model is computed and compared with that one 
obtained with a foamy space modeled by N Schwarzschild wormholes. The comparison leads to the 
conclusion that a foamy spacetime cannot be considered as a collection of N Reissner-Nordstrom 
wormholes but that such a collection can be taken as an excited state of the foam. 
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' I. INTRODUCTION 

o : 

Among the long-standing problems in physics, a fundamental one, from a theoretical point of view, is the absence 

of a quantum theory of gravity. Even if a lot of work has been devoted, especially in the direction of string theory 

i-O ■ and subsequent modifications, i.e. the branes, a complete scheme is still lacking. On the other hand in the traditional 

r ^ , path integral approach to quantum gravity, a typical problem is that the generating functional 
M-H 



J V [g^A expiS g [g^] (1) 



is ill defined because T> [g^ v ] does not represent a measure and moreover there is no convergence. However, in the 
context of a WKB approximation with a fixed background, one can obtain interesting informations. Indeed, if one 
£^ ■ considers a background g^ v , the gravitational field splits into 
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O 
<N 
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where S g 2 ^ [h^] is the action approximated to second order. In the context of the Euclidean path integration, last 



9nu = 9nv + h^, (2) 
where h^ u is a quantum fluctuation around the background field. Then Eq. (|l]) becomes 

J Vg^y exp iS g [g^] ~ exp iS g [g^] J Vh^ exp iSf ] [h^ v ] , (3) 



expression becomes 

J Vg„ v exp {-I g [g^]) ~ exp (-I g [g^]) J Vh^ exp (-if* [V,]) . (4) 



Previous equation can be cast in the form 



r = ^exp(-J d ), (5) 



where T is the decay probability per unit volume and time, while A is the prefactor coming from the saddle point 
evaluation and I c i is the classical part of the action. If a single negative eigenvalue appears in the prefactor A, it means 
that the related bounce shifts the energy of the false ground state ffl. In particular, it is possible to discuss decay 
probabilities from one spacetime to another one H-fl]. For certain classes of gravitational backgrounds, namely the 
static spherically symmetric metrics, it could be interesting the use of other methods based on variational approach. 
In a series of papers, we have used such an approach to show that a model of spacetime foam can be concretely realized 
if one considers a collection of Schwarzschild wormholes whose energy is given by a Casimir energy [p|-|l3| . We recall 
that the Casimir energy involves a subtraction procedure between zero point energies having the same boundary 
conditions. In this paper we would like to apply such methods to another class of static spherically symmetric 
wormholes involving the electromagnetic field. This class is described by the Reissner-Nordstrom (RN) metric which 
is a solution of the Einstein-Maxwell equations with two basic parameters: a mass M and a charge Q. Together with 
the Schwarzschild metric also the RN metric shares the property of being asymptotically flat. This means that the 
Casimir energy computation for such wormholes class will be made with flat space as a reference space. The final 
expression will be compared with that one obtained for the Schwarzschild wormholes showing that the Casimir energy 
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for RN wormholes is always higher than the Casimir energy for the Schwarzschild wormholes. This means that RN 
wormholes cannot be taken as a representation of the ground state of a foamy spacetime. To this purpose we will fix 
our attention to the following quantity 

E RN (M, Q) = E^ (0) + AE™ t (M, Q\ classical + AE™ (M, Q) {1 _ loop , (6) 

representing the total energy computed to one-loop in a RN background. E' lat (0) is the reference space energy, 
which in this case is flat space. AE^ t (M, Q)\ c i assica i is the energy difference between the RN and the flat metrics, 
stored in the boundaries and AE^ t (M,Q)^ 1 _ loop is the quantum correction to the classical term. The analogy 
between AEf^ t (M, Q)u_i oop and the prefactor of Eq.(||) is the key to obtain information about instability. Indeed 
if we discover that the spectrum of the second order differential operator associated with the quantum fluctuations 
of transverse and traceless tensors (TT) admits bound states, then we have instability. In this paper we will assume 
that only one unstable mode appears in the spectrum of TT tensors; this is sufficient to apply Coleman arguments 
about transition from a vacuum to another one. Note that with this assumption we are in the worst situation, namely 
the Schwarzschild and RN wormholes can be compared. In fact without the negative mode a spontaneous transition 
from vacuum to vacuum cannot happen. Therefore in this framework, it is sufficient to compute the stable part 
of the Casimir energy to compare these different pictures. To concretely compute Eq.(Q) we refer to the following 
Hamiltonian with boundary 

H T = H^+ Hox = ( d 3 x (NH+NiH*) + H BS , (7) 
where N is the lapse function, Ni is the shift function and 



H* = -2tf.+W M . 



(8) 



Wm is the energy-momentum tensor contribution of the electromagnetic field, -ffas represents the energy stored into 
the boundary. Since the metric considered has no off-diagonal elements, the Hamiltonian becomes 

H T = Hz+H ds = [ d 3 xNH + H 9s . (9) 



E 



However the physical quantity of interest is the Casimir energy. Therefore we will consider the following expectation 
value 



Hg N - Hl lat 



AEfZ (M, Q) = ± m (10) 



where H2F and H*} at are the total Hamiltonians referred to the RN and flat spacetimes respectively for the volume 
term and * is a trial wave functional of the gaussian form. Note that the flat space is a particular case of 

the RN space with the parameters M = Q = 0. Since the reference space is the same one of the Schwarzschild space, 
it is immediate to recognize that if an instability appears, flat space could decay via a Schwarzschild black hole pair 
creation or via a Reissner-Nordstrom black hole pair creation. Therefore to give indications about the 11 ground state" 
of a foamy spacetime, it is important to establish if 

AEf£ t (M, Q) § AE^ t warzschild (M) (11) 

to one-loop approximation. It is important to remark that it is the vacuum energy difference with the same asymptotic 
reference space that gives the possibility to choose which vacuum is appropriate and not the direct energy difference 
between these possible candidates. The rest of the paper is structured as follows: in section [TJ, we introduce the 



Reissner-Nordstrom metric; in section III, we compute the boundary energy by means of quasilocal energy; in section 
IV , a variational calculation will be set up to compute the Casimir energyand in particular we will restrict our analysis 
to the transverse-traceless tensor metric fluctuations (TT) ; in section [V], we evaluate the spectrum of the operator 
associated with TT tensors (Laplace-Beltrami operator) and we compare the result with what we have obtained in 



case of Schwarzschild wormholes; in section VI, we summarize and conclude. Units in which H = c = k = 1 are used 
throughout the paper. 
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II. THE REISSNER-NORDSTROM METRIC 



The RN line element is 

ds 2 = -/ (r) dt 2 + f (r)" 1 dr 2 + r 2 dft 2 , (12) 



/w .„_!HS + ff,. (13 , 



with 



where Q 2 = G (Q 2 + Q 2 m ); Q e and Q m are the electric and magnetic charge respectively. When the electric charge 
is considered the electromagnetic potential assumes the form A a = (Q e /r, 0, 0, 0) and the electromagnetic tensor 
F a p — d a Ap — dpA a is Fqi — —Q e /r 2 . In the case of a pure magnetic field, the form is A a — (0, —Q m sin 9, 0, 0) and 
the electromagnetic tensor becomes F23 = Q m &in0. Therefore, although the gravitational potential / (r) assumes 
the same form, the gravitational perturbation contribute in a different way. When Q = the metric describes the 
Schwarzschild metric. When Q — M — 0, the metric is flat. For Q ^ 0, we can distinguish three different cases: 

a) MG > Q. In this case the gravitational potential / (r) admits two real distinct solutions located at 




MG+ ^{MGf~Q 2 
MG - (MG) 2 - Q 2 ' 

with / (r) > for r > r + and < r < r_ . In the wormhole language, we will say that r_ is the inner throat 
and r + is the outer throat. In the horizon language r_ is a Cauchy horizon and r + is an event horizon. For 
each root there is a surface gravity defined by 

«± = lim JlflooMI. ( 15 ) 

whose values are 

«+ = (r+-r_)/2r 2 

k_ = (r_ - r+) /2rt [W) 

and for each surface gravity there exists a bifurcation surface associated to a wormhole throat. The Hawking 
temperature associated with the surface gravity of the event horizon is 

T„ = (17) 



b) MG = Q. This is the extreme case. The gravitational potential / (r) admits two real coincident solutions 

2 

located at r + = r_ = r e = MG and its form is / (r) = (1 — MG/r) . Here we discover that k+ = k_ = and 
T H = 0. 

c) MG < Q. In this case the gravitational potential / (r) admits two complex conjugate solutions located at 

r +ti - MG + 1 J Q 2 - (MG) 2 
r_, 4 = MG -i\jQ 2 - (MG) 2 

respectively. 

Cases a) and b) imply Q = when M = 0. We will restrict our attention on case a) only. 
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III. QUASILOCAL ENERGY 



Quasilocal energy is denned as the value of the Hamiltonian that generates unit time translations orthogonal to the 
two-dimensional boundaryQ 

AE™ t (M, Q) ]classical = ^ d 2 x^ (k - k°) , (19) 

where | N\ = 1 at S 2 and k° is the trace of the extrinsic curvature corresponding to the reference space, which in this 
case is flat space. The radial coordinate x continuous on M. is defined by 

dr 

dx = ± (20) 
y/l ^ + # 

which in its integrated form becomes 

/ . It- Md 4- <r . I Ur\ \ \ 

(21) 



:(r) = ±[rVm + MG^[ r - MG + r ^ 

^J(MG) 2 - Q 2 



with x (r) = when r = r± . The surfaces located at r + and r_ are bifurcation surfaces denoted by S+ and 5° , 
respectively. 



r=0 

VWV^ 




r\ / IV \ /r 

VWVvV 

r=0 

FIG. 1. Penrose diagram for the Schwarzschild spacetime. 



In analogy with the Schwarzschild case, see Fig.0, we will restrict our analysis to the regions I and II of FigJ|, 
corresponding to the bifurcation surface S\. The constant time hypersurface we will consider will be denoted by 

S + = £++ U E-| , where the plus sign of Eq.(pb]) is relative to while the minus sign is related to £_| In X++ 

the evaluation of AEp^ t (M, Q)\ c i assica i can be obtained as follows: first we consider the static Einstein- Rosen bridge 



associated to the RN space |17|Jl 



1 See Appendix^] for the explicit derivation of the Hamiltonian in presence of a bifurcation surface. 



4 



FIG. 2. Penrose diagram for the Reissner-Nordstrom spacetime. 

ds 2 = -N 2 (r) dt 2 + g xx dx 2 + r 2 (x) dfl 2 , (22) 

where N, g xx , and r are functions of x defined by Eq.(|2fj|). Second, we consider the boundary S 2 _ + , located at 
x (r) = x ++ (R), and its associated normal = 6£. The expression of the trace 



k = - 



7K (f")s 



(23) 



gives for the RN space 



,V7i 



r \rn 



r \rn 



2MG Q 2 



1 - 



(24) 



Note that if we make the identification TV 2 = 1 — 2M r G + %■ , the line element ( p^ ) reduces to the RN metric written 
in another form; for our purposes the form of N (r) can be left unspecified. Thus the computation of E + gives 



AEPZt (M,Q), las 



8ttG 



dn 2 



(25) 



J \r=R 



R 

G 



IMG Q 2 
R R 2 



When R > 2MG, (and therefore R > Q), we obtain 

AE™ t (M,Q), la 



' \classical 



l - 



MG 
~R~ 



2i? 2 J 



= M. 



(26) 



(27) 



As expected, since the RN space is A.F., we have obtained that the classical contribution to the energy is exactly the 
Arnowitt-Deser-Misner mass (ADM) |Q. Since the RN metric is A.F., the computation of the classical energy term 
leads to 



E RN (M, Q) = E 



Flat 



AE^at (M,Ql clas 



AE™ at (M,Q) 



\clo 



'A J * 



(28) 
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which can be computed by means of quasilocal energy. Thus in complete analogy with the Schwarzschild case, we 
conclude that flat space cannot decay into RN space because the associated boundary energy (ADM) is different, for 

every value of the mass M included the extreme value. However if we look at the whole hypersurface S + = S ++ US H , 

the total classical energy becomes 



E RN (A/, Q) = AE™ t (M, Q)+ lassical + AE™ t (M, Q)^ assical (29) 
with 

AE™ at (M, Q)+ lassical = ^gJ 2 (k k°) , 



(M, Q\ classical = ~ I d 2 xV^ (k - k°) . (30) 
J s._ 



Here the boundaries S+ + and 5? are located in the two disconnected regions M.++ and respectively with 

coordinate values x = x ++ (x^ ) and the trace of the extrinsic curvature in both regions is 

kRN = f -2r, / r on E ++ 
\ 2r, x /r on L + _ 



Thus one gets 

'\classical 1 _M on g 



A.E^,(.:w.g,,,,^ , = \ M M on J%+_ , (32) 



where for E_ we have used the conventions relative to £_| and S+_. Therefore for every value of the boundary R, 

(provided we take symmetric boundary conditions with respect to the bifurcation surface), we have 

E RN (M, Q) = M + (-M) = 0, (33) 

namely the energy is conserved. As stressed in Ref. Jl^j, since we have a spacetime with a bifurcation surface, 
the quantities AE$£ t (M, Q)~\ c i assica i and AE R f at (M, Q)~[ classica i have the same relative sign, while the total energy 
is given by the sum AE R % t (M, Q)^ classlcal + AE^ t (M , Q)~ lassica ^ The energies associated to the boundaries 
are symmetric and they have the same relative sign while the total energy reflects the orientation reversal of the 
two boundaries. Apparently, Eq.(]30|) seems to violate causality, because of the minus appearance in front of the 
expression. However, as reported in Appendix^, it is such a sign that makes the computation causality preserving. 
Since the total classical energy is conserved we can discuss the existence of an instability. To this aim we refer to the 



variational approach to compute the energy density to one-loop [14-16 2C|,|21| 



IV. ENERGY DENSITY CALCULATION IN SCHRODINGER REPRESENTATION 

In previous section we have fixed our attention to the classical part of Eq. (^|) . In this section, we apply the same 
calculation scheme of Refs. |l4|-|l6| to compute one loop corrections to the classical RN term. Like the Schwarzschild 
case, there appear two classical constraints 

Wll (34) 
and two quantum constraints 

' H * = ° (35) 



2 In Ref. |17| we have a subtraction instead of a sum. This is due to conventions adopted. 



G 



for the Hamiltonian respectively, which are satisfied both by the RN and fiat metric on shell. = is known as 
the Wheeler-De Witt equation (WDW) . Our purpose is the computation of 



ttRN 



flat 



<*!*} 



(36) 



where H§ N and H^ lat are the total Hamiltonians referred to the RN and flat spacetimes respectively for the volume 
term || and "J is a wave functional obtained following the usual WKB expansion of the WDW solution. In this context, 
the approximated wave functional will be substituted by a trial wave functional of the gaussian form according to the 
variational approach we shall use to evaluate AEj^ t (M, Q)\i_i oop - To compute such a quantity we will consider on 
£ deviations from the RN metric spatial section of the form, 



+ h 



(37) 



where 



l , il ,ly,l.r l = ( 1 - + dr 2 ; r-.M 1 



(38) 



is the spatial RN background. Thus the expansion of the three-scalar curvature J d 3 Xy/gR^ up to o (ft 2 ) gives 



/ d 3 x^ 



-~hAh + h h Ah H - ~h ij ViVih l j + -WiVifc" - ^R ia h a 3 + hiR^W 



f dWl [z /l ( i * ( ° ) ) /l 



(39) 



where R^ is the three-scalar curvature on-shell. To explicitly make calculations, we need an orthogonal decomposition 
for both TTij and hij to disentangle gauge modes from physical deformations. We define the inner product 



(h, k) := / x TiC'"!,,, (x) k kl (x) d 3 x, 



(40) 



by means of the inverse WDW metric Gijki , to have a metric on the space of deformations, i.e. a quadratic form on 
the tangent space at h, with 



Qijkl = ( g ik g jl + g il g3 k _ 2 gVg kl ). 

The inverse metric is defined on co-tangent space and it assumes the form 

(P, q) ■= / \^ //<''' (x) q kt (x) d 3 x, 

so that 

G« nm G nmkl = \ (Si6f + SfSi) ■ 



(41) 



(42) 



(43) 



Note that in this scheme the "inverse metric" is actually the WDW metric defined on phase space. Now, we have the 
desired decomposition on the tangent space of 3- metric deformations |22j23j : 



1 



hgi 3 + (L£) 



is ■ ^ 



where the operator L maps & into symmetric tracefree tensors 

(L^ij = Vi^ + V 3 -& - | 9y (V • . 
Then the inner product between three-geometries becomes 



(44) 



(45) 
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(h, h) := f x gG ijkl hij (x) h kl (x) d 3 x = 



V9 



tf + iLZfiLQy + h^hjj 



With the orthogonal decomposition in hand we can define the trial wave functional 



* [hij (it)] = Afexp 



Trace 

y 



(46) 



(47) 



where M is a normalization factor. Since we are only interested in the perturbations of the physical degrees of freedom, 
we will only fix our attention on the TT (traceless and transverseless) tensor sector, therefore reducing the previous 
form into^l 



(48) 



This restriction is motivated by the fact that if an instability appears this will be in the physical sector referred to 
TT tensors, namely a nonconformal instability. This choice seems to be corroborated by the action decomposition of 
p4| , where only TT tensors contribute to the partition function^} To calculate the energy density, we need to know 
the action of some basic operators on \P [hij] p0|| . The action of the operator hij on |$) = 5» [hij] is realized by 



while the action of the operator 7Ty- on | , in general, is 



Try (a:) |*) 



Shtj (J?) 



* [hi. 



The inner product is defined by the functional integration 

(*i | *a> = / Phij] *X{hij} * 2 {h k i} 



(49) 



(50) 



(51) 



and by applying previous functional integration rules, we obtain the expression of the one-loop-likc Hamiltonian form 
for TT (traceless and transverseless) deformations []l4|--|l6[ 



1 



M 



ijkl 



i 



(IGttG) (x, x) ijM + — (A 2 )« K± (x, x) 



' iakl 



The propagator K {x,x)iaki comes from a functional integration and it can be represented as 



ia ki 



h 2A(r) ' 



(52) 



(53) 



where h^ a (~x*) are the eigenfunctions of A 2 - t denotes a complete set of indices and A (r) are a set of variational 
parameters to be determined by the minimization of Eq. (p2[) . 



3 In this paper we have defined the trial wave functional without a Planck length constant factor in the exponent. This choice 
alters momentarily the physical dimensions of the problem which are reestablished after the variational procedure. 
4 See also pBl for another point of view. 
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V. THE REISSNER-NORDSTROM METRIC SPIN 2 OPERATOR AND THE EVALUATION OF THE 

ENERGY DENSITY 



To evaluate the energy density, we are led to study the following eigenvalue equation 

(A 2 )]hi = E 2 h) (54) 
where (A2)" is the Spin- two operator for the RN metric defined by 

(A 2 ); := -AS- + 2Vf. (55) 
A is the curved Laplacian (Laplace-Beltrami operator) on a RN background and V? is defined as 

17- := - /';••. (56) 
where R°j is the mixed Ricci tensor whose components are: 

J.^ + **^) (57) 



and Pj is a mixed tensor coming from the electromagnetic energy-momentum tensor expanded to second order in /i^. 
We will follow Regge and Wheeler in analyzing the equation as modes of definite frequency, angular momentum and 
parity |pq| . The quantum number corresponding to the projection of the angular momentum on the z-axis will be set 
to zero. This choice will not alter the contribution to the total energy since we are dealing with a spherical symmetric 
problem. In this case, Regge- Wheeler decomposition shows that the even-parity three-dimensional perturbation is 



diag 



H(r)[l-^^- + %) , r z K(r),r' sin' $K(r) 



Y m (0,<t>). (58) 



In this representation H (r) and K (r) behave as they were scalar fields and the Laplacian restricted to S is 

/ 2MG Q 2 \ d 2 ( 2r-3MG Q 2 \ d 1(1 + 1) 
A; = 1 — ~jT 1 ~Tj2 \ la + T X 15—- ( 59 ) 



The mixed tensor in Eq. (|36|) assumes a different form according to whether we are dealing with the magnetic or 
electric charge. 

A. The electric charge contribution 

In this case the mixed tensor becomes 

P/ := [T aP u a uP) 5% - (F 01 ) 2 515] (60) 



and for a generic value of the angular momentum L, the system ( |54j ) becomes 

(61) 



-A, - ^ + W£) H t (r) = Elf H Hi (r) 



-A l+ 2 -^- Ki (r) = Ef K Ki (r) 
where Ef H and Ef K are the eigenvalues for the Hi (r) field and the Ki (r) field respectively. Defining reduced fields 



5 See Appendix ^. 
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Hi (r) = tf,(r) = ^, (62) 

r r 

and passing to the proper geodesic distance from the throat of the bridge defined by Eq. ( p0[ ) , the system ( |6l[ ) becomes^] 

-a^&l (*) + ^ + (as) (a;) = 

(63) 

-£h (x)+Vr (x)h {x) = Efk{x) 1 



where 



^ w _£m T «« ± «L + ML. (64) 

(x) r J (x) r 4 (x) r 4 (or) 



When r r + x (r) ~ r, then we can approximate the potential with 



^(^^[£+11 (65) 



and the solution for the system (pj) is 



/ipi (a;) = fc p ; (x) = y ~ (px) j/ (pa;) , (66) 
where ji (px) is the spherical Bessel function. If we consider flat space, i.e. M = Q = 0, the system ( |63| ) becomes 

(-£ s + l -^)h l (r)=Efh l (r) 

(67) 

+ W = ^ (r) 



and the solution is 



h p i (r) = k p i (r) = J ^ (pr) jj (pr) . (68) 



On the other hand when r — > r > r_( 



/(/ + !) 3MG 6Q 2 3g 2 



x (r) ~ ^2k+ (r - r+) (a .) _^ v y T — g- ± + = const > ( 69 ) 

where k + has been defined in Eq. (|l6|). However to use the simplicity of Bessel functions for flat and curved space 
when r — > ro > r+, we approximate the potential with 



(x) = ^^ T ^±^ + ^. (70) 



t/t / \ + ^ 3MG 6Q 2 3Q 

X /q /q /g 

Thus close to the wormhole throat we experience a potential barrier (potential hole) whose solution, for system (|6 
is 

| hp_l{ X )= yf^(P_ X )jl{P-X) 

P T is such that P| = £ 2 =F c 2 whereFPI 



(71) 

fcp +( (^)= J^{P+x)ji{P+x) 



The system is invariant in form if we make the minus choice in Eq. (|20[) . 

7 Actually a more standard approach of this problem can be considered by means of phase shifts Si (p). 

8 This choice is also dictated by the necessity of avoiding that different values of the angular momentum enter in the approximate 
potential V* fx) couple like in Refs. §,Ql|. 
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a _ 3MG QQl 3QI 

C= F r 3 4~"~ 4 ■ 

rg 

Note that for functions described by Eq.(|6^) or Eg. (|7l|) we have that 

hi (x) , ki (x) — > when a; (r) — > x (r+) ~ 

Thus the propagator becomes 

' fc P+; (x) k P+l (y) /A+ (P) 



PJ± (x, y) = 1 V 
r fx) r (v) ^— ' 



(72) 



(73) 



(74) 



A± (P) is referred to the potential function V ; (x). This is the more general expression for the propagator. Indeed, 
when one considers flat space or the region far away from the throat r + , it is sufficient the substitution of Pip with p. 
Inserting Eq.(^3) into Eq.(j52|) one gets 



2-7T ^— ' 47T 

p,i i=l 



(16ttG) Ai(p) 



Ef(p,M,Q e ) 



(16ttG) a, Op) 



(75) 



where 



P^(p,M,Q e )=p 2 ~3MG + ^ 



E 2 (p,M,Q e )=p 2 + 



2 i 3MG 3Qf 



(76) 



Ai (p) are variational parameters corresponding to the eigenvalues for a (graviton) spin-two particle in an external 
field. By minimizing f j7^ ) with respect to (p) one obtains A, (p) — [E 2 (p, M, Q e )l 2 / (16ttG) and 



P ,i 



Ef(p,M,Q e ) + JEi(p,M,Q e ) 



(77) 



with 



p 2 + ct > 0, 



(78) 



where we have used definition (|72j) . We can evaluate the total energy for the electric RN background by replacing the 
sum with an integral leading to 



E(M,Q e ) = ^J o °°dpp 2 ^ 



+ ct + \ p l + c 



where V is the volume localized near the wormhole throat. For flat space we put M = Q = and we get 

V rca 



P(0) 



dpp 6 



(79) 



(80) 



Now, we are in position to compute the difference between (|79j) and (j8Cj). The explicit evaluation of the integrals of 
Eq.(f§) in the U.V. limit, gives 



E (M, Q e ) 



V 1 
4^4 



2 P 4 +P 2 (4 



c 2 ) - -± In 



P 



(81) 



Thus 



AEfZ (M, Q e ) = E (M, Q e ) - P (0) 



4^2 



dpp 



p 2 + c 2 . 



+ 4 



2 VP 2 
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V 

4^ 



2 I 2 2\ 

-p (c + + c_) - 



1G 



In 



A 2 



(82) 



where we have used the approximation p 2 >> c 2 ^ and a cut-off A < m p to keep under control the UV divergence 
and we have introduced an arbitrary scale /i. In particular for the Schwarzschild case, [i has been determined by the 



quantity 3MG/r% 



3 _ „2 
-M 



and in that case the approximated expression for AE (M) is 



AE (M) ^d^ f*? 

2ir 2 16 1 " 2 



3 A 2 



32tt 2 V r o / \3MG 



(83) 



Nevertheless, in this paper we would like to establish if a mechanism of space-time foam formation can arise in 
competition with the foam model created by Schwarzschild wormholes. To this purpose in this case we shall fix as a 
scale the value c 2 M and we define the dimensionless parameter a = (Q e /MG). Thus Eq. (Q) becomes 



AEfZ (M, Q e ) ~ — 2 K 



V 



1 

To 



[l + 5a 4 e ~ 4a 2 e ) x 2 In: 



(84) 



where we have defined 



,MG 



3MG 



ro ' r 3 A 2 



Since MG > Q and since r > MG > Q then < a 2 e < 1. In Fig.||, we show the plot of Eq.(|8 
of the parameter a 2 . 



(85) 



for different values 
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5x10 
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FIG. 

space. 



3. Plot of Electric charge contribution for different values of the parameter a e . For a e = 0, we have the Schwarzschild 



Since only the minimum of AE^ t (M, Q e ) is physically relevant, we see that for a 2 = 10 _1 , the minimum is positive. 
This means that the electric field tends to dominate the gravitational attraction and becomes more and more stronger 
as a 2 approaches the extreme value, i.e. MG = Q e . As expected for small values of a 2 , the behavior for the energy 
gap of the Schwarzschild and RN backgrounds with flat space as reference space are very close. Nevertheless, the 
energy gap of the last one is always higher than the same one computed with the Schwarzschild background. This 
means that even for very small electric charges, the Casimir energy for RN wormholes is always higher than the 
Casimir energy for Schwarzschild wormholes. 



B. The magnetic charge contribution 



For the magnetic case P a is 



T, 



a/3'' 



2(P2 3 ) 2 (6$% + S$5$ 



(86) 
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Since the last term in previous equation mixes the components h\ and h%, we have 



2MG 6Q 2 r 



I 2 + 



1 

1 



hi 



(87) 



By repeating the same steps of section V A , the system (|54|) becomes 



-A, 



4MG 



^\ H (r)=E 2 H H(r) 



where we have diagonalized the operator in (|S7|). Since the diagonalization gives the eigenvalues ±1 with eigenvectors 

^2_j_^3 ^ _ .... . . /i 2 

2 ^ 3 respectively, we have considered only the positive eigenvalue since the negative one has the eigenvector ^ 3 

which vanishes because h\ — h\ in Regge- Wheeler representation. Thus following the steps we have used for the 
electric part, we arrive to 



.^h(x) + V+ (x)h(x)=Efh(x) 
-£,k{x)+V l - (x)k(x)=Efk(x) 



where 



V* (x) = 



1(1 + 1) 3MG 4Q 



( x ) r (x) 3 r (x) 4 r (x) 4 



(89) 



(90) 



In analogy with the electric case, we obtain the approximate solutions of system (|39|), by restricting the analysis to 
the sector where r — ► rg > r + .Then 



h P _l{x)= yJ%(P-X)jt(P-X) 

k P+l {x)= JI(P+x)ji(P+x) 



(91) 



P T is such that Pi = Ef =F <& , where 



3MG AQ 2 m , hQl 



Finally we arrive to 



with 



_ V 2 f°° I 

E (M, Q m , A) = Yl / d PP 2 V E ? fa M ' g*0 



EUp,M ) Q m )=p*-mG + !% 



E$(p,M,Q r , 



3MG , 9<3„ 



with the usual condition 



p 2 + 4 > 0, 

where we have used definition (p3). For the magnetic RN background we get 



E(M,Q m ) = -^ J dpp 2 (jp 2 + d 



The zero point energy for flat space is given by Eq. 
like the electric case, we find 



p* + d\ 

then the Casimir energy is given by Eq. 



(92) 
(93) 

(94) 

(95) 

(96) 
Proceeding 



13 



V 



AE™(M,Q m )^—A 



a~ I x 2 In x 



(97) 



where we have defined = a 2 AIG/ro and x is the scale variable we have used for the electric case and < a^i < 1. 
In Fig.|[ we show the plot of Eq. (|?]) for different values of the parameter a^ n 



5x10 



0.2 OA 0.6 0.8 1 



FIG. 4. Plot of Magnetic charge contribution for different values of the parameter a^. For — 0, we have the Schwarzschild 
space. 



Even for the magnetic case, we find a behavior of the energy gap similar to the electric one. This means that 
even in this case, the Casimir energy for RN wormholes is always higher than the Casimir energy for Schwarzschild 
wormholes. 



VI. SUMMARY AND CONCLUSIONS 



In this paper we have considered the possibility of forming a foamy spacctime using a collection of N coherent RN 
wormholes. By applying the same methods used for the Schwarzschild wormholes, we have found that in case of a 
electric or magnetic potential the Casimir energy is always higher than the one computed for the Schwarzschild case, 
i.e. 

AEf& (M, Q) > AE^ a h t warzschUd (M) (98) 

. To this purpose to correctly compare such energies the same reference scale has been considered. We would like 
to remark that for very small electromagnetic field energy contribution, it is the total energy difference between the 
RN spacetime and flat space which is negative. The single energy contribution of flat space and RN space is strictly 
positive. Note that this result has been obtained when the number of wormholes is equal to one. Indeed the complete 
conclusion can be supported when one proofs: 

a) the existence of an unstable mode; 

b) a boundary reduction mechanism coming into play in such a way to stabilize the system. 

Nevertheless, the purpose of the paper is to proof that the Casimir energy for RN wormholes is always higher than 
the Schwarzschild Casimir energy. Since this is the case, it is sufficient to assume the validity of point a) and b) to 

2 

conclude^. Note that the result is obtained for every value of a\ and c? m s.t. < f < 1. This seems to suggest that 

a space-time foam formation realized by RN wormholes is suppressed when compared with the foamy space formed 
by Schwarzschild wormholes. On the other hand one can think to the collection of N RN wormholes as an excited 
state with respect to the collection of N Schwarzschild wormholes leading to the conclusion that such a collection can 
be considered as a good candidate for a possible ground state of a quantum theory of the gravitational field, when 
compared to a superposition of large N RN wormholes. 



9 For the point a), the existence of an instability has been proved in Ref. [£7J only for RN solutions with magnetic charge. As 
far as we know the case with electric charge seems to be stable. Nevertheless, it is interesting to note that the AdS4-RN shows 
such an instability M 
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APPENDIX A: KRUSKAL-SZEKERES COORDINATES FOR RN SPACETIME 

We have defined the RN line element in Eq.(|lJ). To introduce the Kruskal-Szekeres p9|-|3T| type coordinates we 
consider the following transformation 



, , , , . V» 

ds 



1 _ ?M£ + 91 \ [ d e - dr* 2 } + rW 



l-^^ + 9l) d vdu + r 2 (u,v)dn 2 , (Al) 



where v — t + r* is the ingoing radial null coordinate and u = t — r* is the outgoing radial null coordinate. The 
"tortoise coordinate" r* is defined by 

, * r 2 dr , , 

dr* = A2 

[r — r_ j (r — r + j 

=»r*=r+ J-ln |r ~ r+l + J-ln'^^'+c, (A3) 
2/t + r + 2k- r_ 



where k + and k_ have been defined by Eq. Qlq) . To avoid singularities we can define Kruskal-Szekeres type coordinates 

V + = expn + v U + = — cxp— k + u. (A4) 

These coordinates do not cover r < r_ because of the coordinate singularity at r = r_ (and U + V + is complex for 
r < r_), but r — and a similar four regions are covered by the (U~, V~) Kruskal-Szekeres- type coordinates to this 
case. Thus, let us define 

V~ = expK_t> U~ = — cxp— (A5) 

For the + sign we have 

U + V + = - exp (k+ (v -u)) = - exp {2n+r*) = - cxp (2n + r) f - " (A6) 

and the respective line element is 

ds\ = - ^ CXP ( ~ 2K+r) {^—] ~ 1 dU + dV + +r 2 (U + ,V+)dQ 2 , (A7) 
k + r \r — r_ / 

while for the — sign we have 

U~V~ = - exp (k- (v -u)) = - exp (2^_r*) = - cxp (2n^r) ( T ~ - T ) ( — -) + (A8) 



and the associated line element is 

-l 

) 

r + — r 

The conformal Penrose diagram of the RN space is shown in Fig.|[ 



dsl = _ r_±^ exp(-2._r) . . , + ^ ^ ^ ^ (Ag) 
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APPENDIX B: THE HAMILTONIAN CONTRIBUTION OF THE ELECTROMAGNETIC FIELD 



The form of the hamiltonian for the electromagnetic field can be obtained with the same method used for the pure 
gravitational field. Let u a a normalized time-like vector field s.t. u a u a = — 1. The form of TLm comes from the 
Einstein-Maxwell action and it can be written as 



(Bl) 



where 



T, 



a/3 



-in 



(B2) 



and F Q7 = d QL A 1 — d 1 A a . A a is the electromagnetic potential which, in the case of a pure electric field assumes 
the form A a = (Q e /r, 0,0, 0) while in the case of pure magnetic field, the form is A a = (0, — Q m sin#, 0, 0). Q e and 
Q m are the electric and magnetic charge respectively. Both of them contribute in the same way to the gravitational 
potential of Eq. (|l3|). Since we are interested in electric type R.N. metrics, the on-shell contribution of T a pu a u^ is 



i J_Ql 

8tt r 4 



(B3) 



Nevertheless, this is not the complete contribution when gravitational fluctuations come into play. Indeed up to 
second order in hij, one gets 



<TxH M = 2k / d 6 x^/ 3 g 



(B4) 



where T^u a u f3 = ± (F 01 ) 2 h\h\. Thus Eq. (|?4|) becomes 



2 hi hi 



d A xU 



M 



2k 
'8tt 



drx^/ 3 g 



\h\h\ (F 01 f + \h\h\ (F 01 ) 2 + h 3 3 h 3 3 (F 01 ) 2 



On the other hand, when we consider the magnetic charge, the on-shell contribution of T a ^u a u^ is 



(B5) 



T af) u a u p 



1 Ql 



However for the magnetic part, one gets 

T aP u a u? a FijF i] = F l0 F kl g lk gi l 



(B6) 



(B7) 



then to second order in one obtains 



d 3 xH M = 



2k, 

167T 



d 3 x\Pg 



4 



-h\h\ ( F 23 f- lh%hl {F 23 f~-hlhl (1 . 



,2_ 7, 
4' 



(B8) 



APPENDIX C: ACTION AND HAMILTONIAN 

Here we follow Ref. @ to extract the form of the Hamiltonian in presence of a bifurcation surface with boundaries. 
We will consider for simplicity the case of a single horizon, like the Schwarzschild one, the generalization for the 
Reissner-Nordstrom case is straightforward. We consider two spacelike Cauchy surfaces E' and £". The region lying 
between £' and S" consists of two wedges M + and M_ intersecting at a two-dimensional surface So- The symbol 
£ t ( ± ) denotes the part of S t located in M±. We also denote those parts of £' and S" which are the spacelike 
boundaries of the wedges M± as E' ± and E'^ . The lapse function N is positive (negative) at M+ (M_ ) and equals 
zero at the bifurcation surface. The vector — —Nd^t is future oriented in M + and past oriented in M_. The metric 
is described in Eq. ( |2^ ) and the covariant form of the gravitational action for this foliation with fixed three-geometry 
at the boundaries of M is 



1G 



s = 



2 k 



1 

2k 



A/ 4 



d?z 



'9 



*'(+) 



-9^ 



1 



t" 



d 3 x y/gW K 



1 



I M- K Jt'(-) K JB- 

5ft denotes the four-dimensional scalar curvature, k = 8irG, and 



d 3 x V^7© - S (+) 



-7© - So(-)- 



(CI) 



V=9= \N\Vh 



(C2) 



where a is the determinant of the two-dimensional metric a a b. K is the extrinsic curvature of £ t as a surface embedded 
in M, while k is the extrinsic curvature of the boundaries S embedded on S. Then the trace of the extrinsic curvature 
of the boundaries B as surfaces embedded in M is 



Q = k — npaP , 



(C3) 



where o/ 1 = M a V a / is the acceleration of the timelike normal u^. B is a three-dimensional timelike boundary such 
that B = B + U B- . The spacelike normal n % to the three-dimensional boundaries B is assumed to be outward pointing 
at B + , inward pointing at £U . It is assumed that the integrations are taken over the coordinates x^ which have the 
same orientation as the canonical coordinates (t,x,9,(j)) of the foliation. The negative sign for the integration over M_ 
reflects the fact that the canonical coordinates are left oriented in this region. Besides the volume term, the action S 

contains also boundary terms. The notation f t t/±\ represents an integral over the three-boundary S± at t" minus an 
integral over the three-boundary £± at t' . Under a 3 + 1 spacetime split, the four-dimensional scalar curvature is 



= R (3) + K^K^ - (K) 2 - 2V li {Ku' J 



(C4) 



where R denotes the scalar curvature of the three-dimensional spacelike hypersurface S. By the use of Gauss' theorem, 
the conditions 



u ■ tl\b = 0, u - a = 0, u ■ u = — 1, n ■ n = 1 



(C5) 



and Eqns. (|Cl|)-(|C3j) and (C4), one can rewrite the total action in the form 

S=^~ I d i x</=g~ [R + K^K^ - (K) 2 ] - - f d 3 x^/^k 
2^ Jm + k Jb + 

- 77- / d 4 xV=9 [R + K^K^ - (K) 2 ] - - f d 3 xV^fk . 
The action for stationary solutions finally becomes 

S=- J Hdt , 

with the gravitational Hamiltonian H given by 



H = [ d 3 x (NH + N l U t ) + - [ d 2 xy/aNk - - [ d 2 x^Nk. 
Jt. kJ s , kJ s _ 



(C6) 



(C7) 



(C8) 
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